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Abstract 

We propose a way to remove black hole singularities by using a partic- 
ular nonlinear electrodynamics Lagrangian that has been recently used in 
various astrophysics and cosmological frameworks. In particular, we adapt 
the cosmological analysis discussed in a previous work to the black holes 
physics. Such analysis will be improved by applying the Oppenheimer- 
VolkofF Equation to the black hole case. At the end, fixed the radius of the 
star, the final density depends only on the introduced quintessential density term 
p~f and on the mass. 

It is well known that the concept of black hole has been considered very fasci- 
nating by scientists even before the introduction of Einstein's general relativity 
(see pLj for an historical review). However, an unsolved problem concerning such 
objects is the presence of a space-time singularity in their core. Such a problem 
was present starting by the first historical papers concerning black holes [3J [31 H] 
and was generalized in the famous paper by Penrose [5] . It is a common opinion 
that this problem could be solved when a correct quantum gravity theory will 
be, finally, constructed (see [5] for recent developments). 

In this Letter, we propose a way to remove black hole singularities at a 
classical level. The idea is to use a particular nonlinear electrodynamics La- 
grangian to address this issue. Such Lagrangian has been recently used in 
various analysis in astrophysics, like surface of neutron stars [3 and pulsars [5], 
and also on cosmological [Hj contexts. In particular, we will adapt the cosmo- 
logical analysis in [S] to the black hole case. The analysis will be improved by 
applying the Oppenheimer-VolkofF Equation [TSl to the black hole case. In this 



1 



way, we show that the density singularity of a black hole is removed too. At 
the end, fixed the radius of the star, the final density depends only on p^, a 
quintessential density term which will be introduced in the following, and also 
on the star mass. 

At this point it is very important to provide to the potential readers of the 
Letter some elements that were part of our motivations for using NonLinear 
Electrodynamics (NLED) in the search for physical methods to eliminate the 
black hole singularity in general relativity Those elements provide a good 
guide to understand the formal analysis that we are going to present below. 

As it has been carefully explained in references [71 |H] , the effects arising from 
a NLED become quite important in super-strongly magnetized compact objects, 
such as pulsars, and particular neutron stars. Some examples include the so- 
called magnetars and strange quark magnetars. In particular, NLED modifies 
in a fundamental basis the concept of Gravitational Redshift as compared to 
the well stablished method introduced by standard general relativity [7]. The 
analysis in reference [7] proved that unlike general relativity, where the Gravita- 
tional Redshift is independent of any background magnetic field, when a NLED 
is incorporated into the photon dynamics, an effective Gravitational Redshift 
appears, which happens to depend decidedly on the magnetic field pervading the 
pulsar. An analogous result has also been obtained in reference [S] for magne- 
tars and strange quark magnetars. The resulting Gravitational Redshift tends 
to infinity as the magnetic field grows larger [71 H] , as opposed to the predictions 
of standard general relativity. 

What it is important for our goal in this Letter is that the Gravitational 
Redshift of neutron stars is connected to the mass-radius relation of the object, 
see references [Hil]. Thus, NLED effects turn out to be important as regard to 
the mass-radius relation, and one can also reasonably expect important effects 
in the case of black holes, where the mass-radius ratio is even more important 
than for a neutron star. Then, from a physical point of view, the formal analysis 
presented in the Letter displays a correct procedure to estimate the crucial 
physical properties stemming from NLED effects in presence of super strong 
magnetic fields. In fact, the formal discussion developed in the Letter shows that 
the quintessential density term permits to construct a model of star supported 
against self-gravity entirely by radiation pressure. 

In this sense, the approach in our Letter is similar to the one introduced 
in the recent papers of references [201 [28] , where it has been shown that it is 
possible to have radiation pressure supported stars at arbitrary value of high 
mass. The result is that trapped surfaces as defined in general relativity are 
not formed during a gravitational collapse, and hence the singularity theorem 
on black holes, as proposed in reference [S], cannot be applied. 

It is well known that the conditions concerning the early era of the universe, 
when very high values of curvature, temperature and density were present [T], 
and where matter should be identified with a primordial plasma [H [21 [IHl E] , 
are similar to the conditions concerning black holes physics. This is exactly the 
motivation because the singularity theorem on black holes has been generalized 
to the Universe [H [TOl [11]. In the literature there are various cases where a 
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particular analysis on black holes has been applied to the Universe [1] [TUl [H] 
and vice versa pi [TOl pTj [T4| . 

We will re-discuss the simple case of an homogeneous and isotropic sphere 
that was analysed in the historical paper of Oppenheimer and Snyder [12) and, 
in a different approach, by Beckerdoff and Misner [T^ . See also [T] for a review 
of the issue. In such a case, the well known Robertson- Walker line-element can 
be used, which represents comoving hyper-spherical coordinates for the interior 
of the star [I] , and in terms of conformal time 77 it reads 

ds^ = a{rf}{-dif + dx^ + sin^ xidO'^ + sin^ edip^), (1) 

where a{rf) is the scale factor of a conformal space-time. Using sin^ x ^-re 
choosing the case of positive curvature, which is the only one of interest because 
it corresponds to a gas sphere whose dynamics begins at rest with a finite radius 

In [U [T2I [13] the simplest model of a "star of dust" has been discussed, i.e. 
the case of zero pressure. In this case the stress-energy tensor is 

T = pu®u, (2) 

where p is the density of the star and u the 4- vector velocity of the matter. 
Thus, working with SttG = 1, c = 1 and h—1 (natural units), following [9], the 
Einstein field equations give only one meaningful relation 

(^)+a=3a, (3) 
which admits the familiar cycloidal solution 

a= |I(l + cos77), (4) 
where oq is a constant, that is singular for r) = tt. 

Now, let us realize our modified model. We re-introduce a more general 
stress-energy tensor, i.e. the one regarding a relativistic perfect fluid [HE]: 

T = pu®u — pg, (5) 

where p is the pressure and g is the metric. In this case, by restoring the 
ordinary time with the substitution [1] [S] 

dt = adrj, (6) 

the Einstein field equations give the relation 

(a? = f «^ - 1, (7) 
where dot represents time-derivative. 

We will use the non-linear electrodynamics Lagrangian studied in [9]. It 
describes the Heisenberg-Euler NLED 
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C^ = -\f + ciF^+C2G^, (8) 

where G = ^rjaPf^vF"^ Ff^" , F = F^^F^'' is the electromagnetic scalar and 
ci and C2 are constants [5]. 

Thus, we also use the equation of state 

P=\p- Pf^ (9) 
where (see Eq.(25) of Ref. 0) 

= ycii?4. (10) 

B is the strength of the magnetic field associated to F. In [5] the equation of 
state Eq. (jO)) has been obtained by performing an averaging on electric and 
magnetic fields. The presence of the quintessential density term pj will permit 
to violate the Null-Energy-Condition of the Penrose's Theorem |5j. 

The equation of state ([9]) is no longer given by the Maxwellian value. Thus, 
by inserting Eq. © and Eq. ([TUl) into Eq. one obtains [5] 



which can be solved as obtained by [9] 

This expression is not singular for values of ci > 0. In fact, by using elliptic 
functions of the first and of second kind, one gets a parabolic trend near a 
minimum value in t = tf for a(t) , where i/ is the time which corresponds to 
the conformal time i] = tt. 

In concrete terms, by calling l,m,n the solutions of the equation SciBq — 
Blx + 'ix^ = 0, it is [9] 



t = [-(m-O^ A(arcsm \/ — — ^, ^ j—)+n{m~l) ^B(arcsm — — ^, 

(13) 

A{x,y)^ I dz[(l-z2)-i(l-yV)-i] (14) 
is the elliptic function of the first kind and 



where 

/•sin X 



sm X 



B{x,y)^ dz[{{l-z')-')-H{l-y'z')-')-^] (15) 

^0 

is the elliptic function of the second kind. 
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Then, by defining the minimum value of the scale factor like 



«/ = <t = tf), (16) 

and recalling that the Schwarzschild radial coordinate, in the case of the 
black hole geometry ([H, is [T] 

r = asinxo, (17) 
where xo is radius of the surface in the coordinates (H]), one gets rj = 
Uf sinxo > 2Af if Bq has an high strength, where M is the black hole mass and 
2M the gravitational radius in natural units 1 . Thus, we find that the mass of 
the star generates a curved space-time without event horizons. 

Now, let us consider the famous Oppenheimer - Volkoff Equation [TS] (r the 
radial coordinate) 



dp _ {M + 4TTr^p){p + p) 
dr ^ r(r - 2M) ' 

Using eq. (O one gets 



Pfi = P04 + r{r-2M] dr Pf5 = P05 " ,(,_2M) 



dr 



the integral (|T9| can be computed by separating variables. 
One then gets 

P/ - Poi ^ ^ + P02 + P.9M + 
+{p„Vl6^[ro(ro - r/) + ^{rl rj) + In ^1}-^+ 
+P04 exp607r[ro(ro - r/) + i(r2 - rj) + 4AP In ^;^] + 



(18) 



where r/ = r(t/), ro = r(0), p/ = p{tf), po = p(0). 
By defining 

Pf = Pfl + Pf2 + Pf3 + Pf4 + Pf5 Pfl = POl - Q r{rf2M) dr 

Pf2 = P02 + r{r-2M) dr Pf3 = P03 " ,(,^2M) ^'^ (20) 



(21) 



+P05 + 367rp2 [ro(ro - r/) + i(r2 - r^) + In '^^]. 

A reasonable assumption is p — s> for r ^> ro, i.e. at the surface of the star. 
In this case it is poi = P02 = Pf 03 = Po4 ~ P05 — 0, and the final density reduces 
to 
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Pf — p-y'^M In 




rj)+iM^ In 




r/ ro — 2M 



(22) 



Thus, we have shown that the density singularity has been removed too 
and that, fixed the radius of the star, the final density depends only by the 
introduced quintessential density term p^ and by the mass. 

Summarizing, adapting to the black hole cosmological analysis in 

[9] , we have shown that the introduction of the Heisenberg-Euler NLED ([8|) in 
the black hole physics can, in principle, remove the black hole singularity for 
values of ci > in Eq. The analysis has been be improved by applying the 
Oppenheimer - VolkofF Equation [T5] to the black hole case. In this way, the 
density singularity of a black hole has been removed too. We found that at the 
end, fixed the radius of the star, the final density depends only by p-y and by 
the mass of the black hole. 

We are going to further improve the analysis by using the Born-Infeld La- 
grangian in the context of black holes physics in a future work [16] . 

Finally, we take the chance to stress that this result comes to follow in the 
same direction of the papers on this issue pioneered by Mitra [171 HB [H] and 
Robertson, Leiter and Schild [H [23l [H [251 [26l [27] , who have shown that the 
formation of trapped surfaces during the gravitational collapse can be avoided 
if the collapsing star develops a powerful magnetic field of strength of > IQ^^G, 
or equivalently, an intrinsic magnetic moment on the order of lO^^G — cm'^, in 
the case of galactic black hole candidates. 

For a sake of completeness [29], we recall previous works attempting to 
remove black holes singularities using nonlinear electrodynamics by Ayon-Beato 
and Garcia in [SO] |3ll|32l|33l |34] . 
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